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METHOD FOR DETERMINING NOISE IN RADIOGRAPHY 
CROSS-REFERENCE TO RELATED APPLICATIONS 

10001] This application claims the benefit of a priority under 35 USC 1 19 (a)- 

(d) to French Patent Application No. 02 13566 filed October 29, 2002, the entire 
contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

|0002] The invention and embodiments thereof is directed to a method for 

determining noise in radiography and, in particular in fluoroscopic noise. The field of 
the invention is that of imaging, and more precisely that of the reduction of 
fluoroscopic noise in images. The field of the invention is more particularly for the 
reduction of noise in images acquired in time-based sequences in radiography. 

|0003] In the prior art, it is known that there is a fluoroscopic noise in 

fluoroscopy and that the mean standard deviation of this fluoroscopic noise is 
proportional to the square root of the number of photons reaching the detector. 

]0004] In the prior art, the presence of this noise is known but is not 

specifically dealt with except on the basis of the experience of an observer using the 
images in order to interpret them. There is then a problem of interpretation, but also a 
problem of exposure to an object to be imaged, such as a patient. It may become 
necessary to increase the radiation dose in order to improve the signal-to-noise ratio 
of the image, i.e., to increase the proportion of relevant information in the image. 
This increase in radiation is not free of consequences for the object, particularly a 
patient, which then has a risk of being destroyed or damaged or subject to greater than 
recommended exposure to radiation. 
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BRIEF DESCRIPTION OF THE INVENTION 

(00051 In an embodiment of the invention, the problem of fluoroscopic noise 

is resolved by providing a model of the noise. When the noise has been modeled, it is 
possible to eliminate it from the image and thus improve the readability of the image 
and, if necessary, reduce the radiation dose required to obtain a readable image. 

|0006] In an embodiment of the invention a method for determining noise in 

radiography comprises: 

10007] acquiring at least two images, i-1 and i, of a same zone; 

]0008] coding the acquired images into digital images that can be identified 

with matrices having horizontal by vertical dimensions equal to N x M, each digital 
image being then formed by N.M dots, each dot of an image i being identifiable by its 
coordinates 0 < x < N, and 0 < y < M, this dot then being referred to as a dot Pi(x,y), 
each dot Pi(x,y) then having a corresponding value v which is the result of the 
acquisition of the image, the value v having a dynamic range from Vmin to Vmax; 

(0009] dividing the dynamic range Vmax - Vmin into sub-groups defined by a 

lower limit Bi and an upper limit Bs, the sub-groups having a null intersection, the 
joining of the sub-groups covering the dynamic range Vmax - Vmin, a dot of an 
image i then belonging to a given sub-group when Bi <= Pi(x,y) <Bs, where Pi(x,y) is 
the gray level of the pixel of the image i with coordinates (x,y); 

]0010] computing, for at least one sub-group SG, of the mean standard 

deviation o of the values Pi(x,y) - Pi-l(x,y); 

(00111 discriminating the values Pi(x,y) of SG to keep only those values such 

that the criterion C : Pi(x,y) - Pi-l(x,y) < u(Pi(x,y) - Pi-l(x,y)) + k.a, is met and thus, 
a sub-group SG' is obtained, where u. is a mean value; 
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[0012] applying the same processing operations to the sub-group SG f as to the 

sub-group SG by iteration until a sub-group SG", corresponding to an end-of-iteration 
criterion, is obtained; 

[0013] performing iterative processing operations on all the sub-groups 

defined in the dynamic range Vmax - Vmin and thus, for each sub-group, a mean 
standard deviation, associated with an x-axis value v = (Bi + Bs)/2, is obtained; and 

[0014] performing an operation of regression on the dots obtained at the 

previous step to determine the coefficients a, P and y of the noise function: a(v) = 
a.Vv + p.v + y defining the noise for a given value v. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] The invention will be understood more clearly from the following 

description and from the accompanying figures, which are given purely by way of an 
embodiment and in no way restrict the scope of the invention. Of these figures: 

[0016] Figure 1 illustrates the method according an embodiment of the 

invention; 

[0017] Figure 2 illustrates an image sequence; and 

[0018] Figure 3 is a graphic representation of the pairs (v, a). 

DETAILED DESCRIPTION OF THE INVENTION 

[0019] In an embodiment of the invention, the problem of fluoroscopic noise 

is resolved by filtering the image obtained so as to improve its quality. Thus, the 
fluoroscopic noise is eliminated after it has been determined/modeled. This noise is 
totally determined, hence modeled, by its mean standard deviation. This mean 
standard deviation is itself a function of the square root of the number of photons 
received by a detector. The number of photons itself is related to a gray level in a 
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digital image. Digital images are thus used to obtain the modeling. This modeling is 
done in several steps. In a first step, two digital images of the same zone are acquired. 
In an embodiment of the invention, an image, unless otherwise indicated, is a digital 
image. It will be noted here that the teaching of the invention is valid whatever the 
nature of the sensor, whether digital or analog, used to obtain the images. Each pixel 
or dot of an image is paired with a pixel of the other image by means of its 
coordinates in the image. Each pixel also has a gray level value, or gray level. The 
pixels are grouped together by gray level intervals and thus sub-groups of pixels are 
obtained. For each sub-group of pixels, a discrimination is made as follows: the mean 
u and the mean standard deviation a of Pi(x,y) - Pi-l(x,y) are computed, where Pi(x, 
y) is the gray level of the pixel of the image i with coordinates (x, y). Then, in a sub- 
group, only the values P(x,y) are kept such that Pi(x,y) - Pi-l(x,y) < u + k. a. This 
discrimination is repeated iteratively on the result of the preceding discrimination. 
This discrimination greatly reduces blur. To obtain even more reliable sub-groups, 
the method eliminates those that, at the end of the discrimination, are not centered, 
namely those whose mean is greater than 1.5 times the mean standard deviation. 
Then there is knowledge of a collection of pairs v, a) where v is a gray level. These 
pairs enable an operation of regression leading to parameters a, P and y such that a(v) 
= a.Vv + p.v + y, where a(v) is the modeling of the fluoroscopic noise. This 
regression is made robust by iterating it after weighting the o values of the pairs (v, a) 
drawing the curve upwards so as to obtain a curve that passes above the majority of 
the dots (v, a). 

[0020] Figure 1 shows a preliminary image acquisition step 101. The 

acquisition is performed by a scanner type radiography apparatus used to acquire 
images of an object, such as internal organs of a living organism, particularly the 
human body. The technologies of these apparatuses are based on the emission of 
radiation received on or at a sensor/detector after it has gone through the object. The 
radiation received by the sensor/detector therefore depends on the object crossed. The 
sensor/detector produces a digital image that can be viewed on a screen, printed, 
and/or processed by a computer, or any apparatus comprising processing circuits. The 
method is recorded in a memory in the form of instruction codes and implemented by 
a microprocessor, in order to improve the interpretation of an image and/or reduce the 
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dose of radiation used to carry out an examination. This memory and this 
microprocessor are preferably contained in the radiography apparatus. In practice, 
this memory and this microprocessor are connected to the radiography apparatus, 
either by an internal connector or by a connector external to the radiography 
apparatus. 

[0021] Classically, an image has 1000x1000 resolution and therefore has one 

million pixels, each pixel being identified by its coordinates in the image and each 
pixel being associated with a value of gray level, intensity or another numerical value. 
In the present example, it shall be assumed that this is a gray level obtained after 
fluoroscopic exposure, namely exposure to weak radiation. An image can be likened 
to a table having N rows by M columns. N and M are not systematically equal to 
1000; their values depend on the spatial resolution of the sensor/detector used. Their 
values do not affect the efficiency of the method. The mode of operation remains 
valid for non-fluoroscopic exposure. 

[0022] In fluoroscopy, the mean standard deviation of the fluoroscopic noise 

is proportional to the square root of the number of photons reaching the detector. 
However, the gray level is proportional to the quantity of photons received. This 
enables working on the gray levels. 

[0023] Fluoroscopic noise is the resultant of the quantum noise and the 

response of the image acquisition system, namely of the detector. The method 
described deals with the totality formed by the quantum noise and the apparatus. 
Technically speaking, the method is not specific to quantum noise or to an apparatus, 
but it can be applied to any other noise that is not spatially correlated. 

[0024] For most of the examinations, the images are acquired in temporal 

sequences, namely one after the other. This makes it possible to observe an object, 
such as an organ during a cycle, for example, the heart cycle.. In an examination of 
this kind, there are therefore several successive images, or shots, of a same region of 
the object. The successive images have the same resolution and two pixels with the 
same coordinates correspond to a same region of the object, unless the object has 
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shifted between the «. shot, Sttch shifts may occur when i, is songh, to obtain an 
image of the heart for example. 

,00251 Consider that the radiography apparatus has acquired two images i-1 

and i of a same region of an object. The images have an N-pixel horizontal resolution 
and an M-pixel vertical resolution, as shown in Figure 2. A pixel with coordmates (x, 
y) of an image i has a gray level equal to Pi(x, y). 

(0026] Step 101 proceeds to step 102 for sub-sampling images coming from 

the radiography apparatus. Each pixel of an image has an associated value 
corresponding to a gray level. This value has a certain dynamic range. In other 
words, each pixel has an associated value ranging from Vmin to Vmax. 
Conventionally Vmin is equal to 0 and Vmax to 2*1* In practice, Vmin and Vmax 
depend on the sensor used and its calibration. When an image is sub-sampled, the 
interval [Vmin, Vmax] is divided into intervals separated in such a way that the 
joining of the sub-intervals covers [Vmin, Vmax]. Each sub-interval [Br, Bs] 

included in [Vmin, Vmax] has a corresponding sub-group of pixels such that Br <= 

Pi(x,y)<Bs. 

|0027] In practice, the sub-intervals have a dynamic range of 20. That is, the 

interval [Vmin, Vmax] is sub-divided into sub-intervals with a width of 20. The sub- 
sampling then comprises distributing the dots of an image in the sub-groups 
corresponding to the sub-intervals. It is possible that there are empty sub-groups. A 
sub-group may also be called a sample. 

[0 028] In practice, the sub-sampling is done as soon as a decision is taken on 

the sub-division into sub-intervals and there is no physical distribute. The 
distribution is done by the simple reading of the gray level of a pixel. The distnbutron 
into sub-groups is done by reading the image i. 

,0029] Step 102 proceeds to step 103 for computing a mean standard 

deviation for a sample. Step 103 considers the dots of a sub-group SG of the image . 
in correspondence with dots of the image i-1. Dots are in correspondence when they 
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have the same coordinates. A sub-group in the image i determines a sub-group 
comprising the same dots in the image i-1. In step 103 a computation is therefore 
made of the mean standard deviation a of the values (Pi(x, y) - Pi-l(x, y)) with (x, y) 
belonging to SG. 

[0030] Step 103 proceeds to step 104 for discriminating the sub-group SG. 

This discrimination comprises the elimination of the aberrant dots resulting for 
example from a blur of the object to be imaged. The discrimination is done according 
to the following criterion: (x, y) of SG belongs to SG' if and only if: Pi(x,y) - Pi-l(x,y) 
< p(Pi(x,y) - Pi-l(x,y)) + k.o, where u(Pi(x,y) - Pi-l(x,y)) is the mean of the values 
(Pi(x,y) - Pi-l(x,y)) for (x, y) in SG. SG" is therefore included in SG but it can happen 
that it does not include all the dots of SG. 

[0031] The number k is a parameter of the algorithm, and is preferably equal 

to 2, but any other non-zero number would be suitable. The greater the number k 
chosen, the less discriminatory is the algorithm. A number k greater than 10 is 
therefore also appropriate but the discrimination is then very weak. 

[0032] Step 104 proceeds to step 105 for determining the end of the iteration. 

A first possible criterion for the end of the iteration is that SG' must be identical to 
SG. In this case, the method proceeds to a centering step 106. If not SG is replaced 
by SG' then steps 103 and 104 are resumed. 

[0033] A second possible criterion for the end of iteration is a number of 

loops. In this case, in step 105, SG is replaced by SG' and the operation is resumed at 
step 103. This replacement is done a predefined maximum number of times, for 
example, five times. However, this number is given by way of an example and can be 
parametrized as a function of the quality of the modeling to be obtained. The greater 
this number, the greater the quality but the longer the modeling. 

[0034] In practice, the first and second criteria can be combined through 

logical combination. The criterion then is that the iteration is done so long as the 
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number of iterations is smaller than the maximum number and so long as SG is 
different from SG'. 

[0035] In step 105, if the end-of-iteration criterion is not valid, the method 

proceeds to step 103 after having performed the operation SG = SG'. If not, the 
method proceeds to step 106. 

100361 Step 106 provides an even more robust method. Step 106 eliminates, 

from the processing, the samples, or sub-groups SG, resulting from step 105 such 
that: u(Pi(x,y) - Pi-l(x,y)) > l.o(Pi(x,y) - Pi-l(x,y)), with (x, y) in SG and 1 being 
conventionally equal to 1.5. The eliminating of a sub-group means that the sub-group 
will no longer be taken into account starting with step 108. 

[0037] Step 106 proceeds to step 107 in which it is verified that all the sub- 

groups determined during step 102 had been processed. If this is not the case, steps 
103 to 106 are repeated for all the sub-groups that have not yet been processed. 

[0038] At the end of step 107, a value v and a mean standard deviation o can 

be associated with each non-empty group that was not removed at step 106. The 
mean value v is the mean value of the sub-interval that was used to initially determine 
the sub-group. Given that a sub-interval is determined by a lower limit Bi and an 
upper limit Bs, v is equal to (Bi + Bs)/2. The mean standard deviation o is the last 
mean standard deviation computed for the sub-group. At step 107 there is collection 
of pairs (v, a). 

[0039] Step 107, when there is no longer any sub-group to be processed, 

proceeds to a regression step 108. Figure 3 illustrates that it is possible to graphically 
represent the pairs (v, o), in placing v on the x-axis, v ranging from Vmin to Vmax. 
Starting from the pairs (v, a), an operation of regression is performed so as to 
determine the first coefficients a, p and y modeling a curve a(v) = a.vv + p.v + y 
passing as close as possible to the set of the pairs (v, a). 
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[0040] Step 108 proceeds to step 109 in which an end-of-iteration criterion is 

assessed. A first end-of-iteration criterion is, for example, the presence of dots under 
the curve defined at step 108. If such dots exist, the method proceeds to step 1 10 in 
which the weighting of these dots is modified. Step 110 proceeds to step 108. If such 
dots do not exist, the method proceeds from step 109 to end step 111. Step 108 is 
therefore performed several times on the basis of a collection of pairs whose 
weighting evolves as a function of the position of the pair with respect to the last 
curve determined, i.e., with reference to the last coefficients a, p and y determined at 
step 108. 

[0041] A second end-of-iteration criterion is the use of an iteration counter. In 

this second criterion, steps 108 to 110 are carried out a predetermined number R 
times. In practice, R is equal to 5. However, it is possible to use any value 
whatsoever for R. 

[0042] An embodiment may use the first and second criteria combined by a 

logic operator, making it possible to take into account the most favorable conditions 
for stopping the iteration in terms of processing time. Thus, the iteration is stopped as 
soon as the maximum number of iterations is reached or as soon as there are no longer 
any pairs below the last determined curve. 

[0043] In step 110, a refining of the regression of the step 108 is prepared. 

This refining is done through a weighting of the pairs (v, a) located above the curve 
defined at the step 108. These dots are easily detected. They are dots for which a(v) < 
a. For these pairs, a weighting P strictly smaller than 1 is applied to the value c of the 
pair (v, a) which then becomes (v, P.a). Step 110 then proceeds to step 108 of 
regression with all the pairs (v, o), whether weighted or not, coming form step 110. A 
typical value of P is 0.95, but any value whatsoever included in the interval [0 ... 
0,99] is appropriate. Satisfactory results are obtained with a value P in the interval 
[0,75 ... 1]. 



9 



14XZ125560 



[0044] When the end-of-iteration condition is fulfilled, the method proceeds 

from step 109 to the end step 111. In step 111, the apparatus having performed the 
processing possesses coefficients a, P and y of the function: a(v) = a.Vv + p.v + y 
modeling the fluoroscopic noise during an examination implementing a fluoroscopic 
noise. Having carried out the processing, the apparatus uses these coefficients either 
to determine a filter enabling the production of an image in which fluoroscopic noise 
has been substantially reduced or to parametrize the fluoroscopy apparatus. 

[0045] In the case of a filtering operation, the image after filtering is therefore 

free of fluoroscopic noise. This favors the reading and interpretation of this image by 
a practitioner. Therefore the practitioner then can place greater reliance on the images 
given by the radiography apparatuses. 

[0046] As for the parametrizing of the fluoroscopy apparatus, it is useful 

because, once the coefficients are determined, the quantity of noise present in an 
image and, hence, the signal-to-noise ratio, are known. If this signal-to-noise ratio is 
satisfactory, i.e., if there is a lot of signal and little noise, the dose of radiation emitted 
by the fluoroscopy device is reduced. The examination then becomes less traumatic 
for the patient. On the other hand, if the signal-to-noise ratio is poor, the radiation 
dose is increased so that the examination is relevant. In both cases, a gain in obtained 
because the optimum dose of radiation needed to obtain the most relevant result 
possible has been successfully determined. 

[0047] In practice, both uses of the coefficients a, p and y are implemented. 

An image sequence comprises a variable number of images. The first two images of 
the sequence are used to determine the coefficients a, p and y that will subsequently 
be applied to the processing of all the images of the sequence. 

[0048] In an embodiment of the invention, the coefficients a, p and y are 

computed for each image of the sequence. Each image is then associated with its own 
set of coefficients. 
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[0049] The description uses an image i to define the sub-groups, and the 

image i is paired with the preceding image i-1. In an embodiment of the invention, the 
image i can also be paired with a following image 1+1. 

[0050] The method and embodiments thereof is implemented either in a 

digital image processing station or in a device for the control of a radiography 
apparatus. With the machines commonly used in radiography, the method and 
embodiments thereof give results within about thirty milliseconds for the processing 
of an image having a definition of one million pixels. These performance levels are 
highly satisfactory in a real-time context. An operator does not feel that he/she is 
waiting for the image. 

[0051] In an embodiment of the invention, the end-of-iteration criterion is a 

number of iterations greater than 5. In an embodiment of the invention, R is in the 
interval [3 ... 10]. In an embodiment of the invention, R is greater than 10. 

[0052] Hence, the method and embodiments thereof provide: (1) for the 

reduction of fluoroscopic noise in images; (2) the reduction of fluoroscopic noise in 
digital images coming from a fluoroscopic acquisition method; (3) the robust 
determining of fluoroscopic noise; (4) the real-time determining of fluoroscopic noise; 
and (5) the regulation of the quantity of radiation received by patient during a 
radiography examination. 

[0053] One skilled in the art may make or propose various modifications to 

the structure and/or steps and/or function and/or way and/or result of the disclosed 
embodiments without departing from the scope and extent of protection. 
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